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Abstract Biological networks, such as cellular metabolic 
pathways or networks of cortico cortical connections in 
the brain, are intricately organized, yet remarkably ro- 
bust toward structural damage. Whereas many studies 
have investigated specific aspects of robustness, such as 
molecular mechanisms of repair, this article focuses more 
generally on how local structural features in networks 
may give rise to their global stability. In many networks 
the failure of single connections may be more likely than 
the extinction of entire nodes, yet no analysis of edge 
importance (edge vulnerability) has been provided so far 
for biological networks. We tested several measures for 
identifying vulnerable edges and compared their predic- 
tion performance in biological and artificial networks. 
Among the tested measures, edge frequency in all short- 
est paths of a network yielded a particularly high cor- 
relation with vulnerability, and identified inter-cluster 
connections in biological but not in random and scale- 
free benchmark networks. We discuss different local and 
global network patterns and the edge vulnerability re- 
sulting from them. 
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1 Introduction 

Extensive evidence shows that biological networks are 
remarkably robust against damage of their nodes as well 
as links among the nodes. For example, Parkinson dis- 
ease only becomes apparent after a large proportion of 
pigmented cells in the substantia nigra are eliminated 
IjDamier et al. 1999(1 . and in spinal cord injuries in rats, 
as little as 5% of the remaining cells allow functional re- 
covery {You et al. 2003|1 . Many metabolic networks, as 
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well, were found to be robust against the knockout of 
single genes. This feature is both due to the existence of 
duplicate genes as well as alternative pathways which en- 
sure that a certain metabolite can still be produced using 
the undamaged parts of the network I ] Wagner 2000| ). 

Approaches of network analysis have been used to in- 
vestigate various types of real- world networks, in which 
persons, proteins, brain areas or cities are considered 
nodes; and functional interactions or structural connec- 
tions are represented as edges of the network {Strogatz 2001| ) . 
Many such systems display properties of small-world net- 
works {Watts and Strog atz 1998), with clustered local 
neighborhoods and short characteristic paths (or average 
shortest paths, ASP). Also, some networks possess more 
highly connected nodes, or hubs, than same-size random 
networks, leading to a power-law degree (scale-free) dis- 
tribution of edges per node, where the probability for a 
node possessing k edges follows fc~ 7 ( Barabas i and Albert 1999 ) 
Such scale-free networks are error-tolerant towards ran- 
dom elimination of nodes, but react critically to the tar- 
geted elimination of highly connected nodes {Albert et al. 2 000 ) 
It has been shown that cerebral cortical networks in the 
cat and macaque monkey brain display a similar behav- 
ior (Kaise r^t al. 20 04). While networks may possess fea- 
tures of both small- world and scale- free organization, the 
two topologies are not necessarily identical. 

Whereas previous studies explored the impact of le- 
sioning network nodes (Baraba si and Albert 1999jl . the 
effect of edge elimination in biological networks has not 
yet been investigated. How can edges that are integral 
for the stability and function of a network be identified? 
In some systems, for instance, transportation or informa- 
tion networks, functional measures for the importance of 
edges, such as flow or capacity, are available. For many 
biological networks, however, such measures do not ex- 
ist. In brain connection networks, for example, a projec- 
tion between two regions may have been reported, but its 
structural and functional strength is frequently unknown 
or not reliably specified (Fcllem an and van Essen 1 991). 
and its functional capacity may vary depending on the 
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task QBuchel and Friston 1997|l . Similarly, in biochem- 
ical networks, reaction kinetics are often highly vari- 
able, or generally unknown (Sch uster and Hilgeta g 1994| 
Stclli ng" et al. 2 002 ) . However, the analysis of a network's 
structural organization may already provide useful infor- 
mation on the importance of individual nodes and con- 
nections, by identifying local features and investigating 
their importance for global network structure and func- 
tion. Specifically, we are interested in these questions: 
Which structural patterns can be identified in a net- 
work? Have biological networks specific features, com- 
pared to artificial networks? As examples for biological 
networks we analyzed cortical fiber networks, metabolic 
networks as well as protein interaction networks and, 
for comparison, also an artificial network, the German 
highway system. 

In our analyses we tested the effect of eliminating 
single edges from networks. For some of the studied net- 
works, such as cortical connectivity, structural damage 
might also result in multiple lesions of connections or 
regions. Such lesions can lead to network fragmentation, 
and may be theoretically analyzed within the framework 
of edge or vertex cut sets ( |Sporns 2002| ). A computa- 
tional algorithm to predict the effect of multiple network 
perturbations (or lesions) was presented as part of the 
MSA (Multi-perturbation Shapley value Analysis) ap- 
proach (Kcina n et al. 2 004). In the present article, how- 
ever, we focus on the impact of eliminating single edges, 
to infer network patterns and their vulnerability. 

Two established parameters for characterizing net- 
works are the average shortest path (ASP) and the clus- 
tering coefficient (CC). The ASP is the average number 
of edges that have to be crossed in order to reach one 
node from another. For a network with N nodes, it is 
calculated as the average of all existing shortest paths: 



ASP = d(i,j) with i ^ j and d(i,j) ^ oo, (1) 

where d(i,j) is the number of edges of the shortest path 
between nodes i and j. Infinite distances between uncon- 
nected nodes, which occur after network fragmentation, 
are excluded from the average calculation. We used the 
deviation of the ASP before and after edge elimination 
as a measure of network damage. 

The clustering coefficient of a node v with k v adjacent 
nodes (neighbors) is defined as the number of edges exist- 
ing among the neighbors, divided by the number k% — k v 
of all possible edges among the neighbors fWatts 1 999). 
We use the term clustering coefficient as the average clus- 
tering coefficient for all nodes of a network. Networks 
with an ASP comparable to randomly connected net- 
works, but with much higher clustering coefficient, are 
called small-world networks ( |Watts and Strogatz 19981 ) ■ 
These systems exhibit network clusters, that are regions, 
in which many interconnections exist within a cluster, 
but only few connections run between clusters. Various 
kinds of networks, such as electric power grids and social 



networks, display small- world properties ( Watts and Strogatz 1 998). 
In addition, neural networks of C. elegans and corti- 
cal networks of the cat and macaque were also shown 
to be small-world n etworks (|Watts and Strogatz 1998"| 
|Hilgetag et al. 2000| |Sporns et al. 2000D and to exhibit 
a clustered architecture ( |Hilgetag et al. 2000| ). 



2 Materials and methods 

2.1 Investigated biological and artificial networks 

2.1.1 Brain networks We investigated long-range fiber 
projections in the cat and macaque brain. Nodes were 
brain regions or areas (e.g., VI), and edges were fiber 
connections between them. We considered connectivity 
data for a non-human primate, the Macaque monkey (73 
nodes; 835 edges; density 16%), and the cat (55 nodes; 

891 edges; density 30%) l|Scannell et al. 1995llgcannell et al. 19991 
|Young 1993| ). In both species, the data included con- 
nections between cortical regions, as well as a few sub- 
cortical structures (e.g., the amygdala) and regions of 
entorhinal cortex. These networks possess properties of 
small- world networks ( |Hilgetag et al. 2000| ) , and also show 
similar response to attack as scale- free networks IjKaiser et al. 2 004 ) . 
The clustering coefficient of these networks ranged be- 
tween 40 to 50% (cf. Table P). 

2.1.2 Protein-protein interactions As an example of bio- 
chemical networks, we examined the protein-protein in- 
teractions of the S. cerevisiae yeast proteome ( |Jeong et al. 200l| > . 
The data consisted of 1,846 proteins and 2,203 distinct 
functional relationships among them, forming 4,406 uni- 



directional edges in the network graph (data from http://www.nd.edu/'- 
With a value of 6.8%, the clustering coefficient was con- 
siderably smaller than for brain networks. As a further 
difference, the protein interactions network consisted of 
149 disconnected components. The main component con- 
tained 79% of the proteins, and the remaining compo- 
nents mostly were composed of only one pair of proteins. 
As shown before ( Jeong et al. 2001] ), connections in this 
network are also distributed in a scale- free fashion. How- 
ever, the yeast two-hybrid method yielding protein inter- 
action data Ijlto et al. 20011 IHavin et aT~2 002) produces 
many artifacts (Kitano 2002) which might have influ- 
enced prediction results. 

2.1.3 Metabolic networks Cellular metabolic networks 
of different species were analyzed. Nodes here were metabolic 
substrates, and edges were considered as reactions l|Ravasz et al. 2002)) 
(Data at http://www.nd.edu/~networks/database/). For 
bacteria, we investigated the metabolism of E. coli with 
765 metabolites and 3,904 reactions. For eukaryotes, the 
data for Arabidopsis thaliana (299 metabolites and 1,276 
reactions) and the yeast S. cerevisiae (551 metabolites 
and 2,789 reactions) were examined. 
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Fig. 1 Frequency of edges in the all-pairs shortest paths and 
resulting network damage after elimination, a Cat brain con- 
nectivity and primate (Macaque) brain connectivity b show 
a strong correlation with damage, c Metabolic network of 
S. cerevisiae. d German highway system. Decreases of ASP 
were caused by eliminated cut-edges, leading to a separation 
of the network. 



2.1.4 Transportation network Data for the German high- 
way (Autobahn) system were explored as a comparative 
example of man-made transportation networks. The net- 
work consisted of 1,168 location nodes (that is, highway 
exits) and 2,486 road links between them (Autobahn- 
Informations-System, AIS, from http://www.bast.de). Only 
highways were included in the analysis, discarding smaller 
and local roads ('Bundesstrassen' and 'Landstrassen'). 
Parking and resting locations were also excluded from 
the set of nodes. Furthermore, multiple highway exits 
for the same city were merged to one location represent- 
ing the city as a single node of the network graph. 



2.1.5 Benchmark networks Twenty random networks 
with 73 nodes and comparable density as the Macaque 
network were generated. Moreover, twenty scale- free net- 
works with 73 nodes and equivalent density were grown 
by preferential attachment (Barabas i and Albert 1 999). 
starting with an initial matrix of 10 nodes. 

In addition to the random and scale- free networks, 
which consisted of only one cluster, networks with multi- 
ple clusters were considered. Twenty networks were gen- 
erated with 72 nodes in order to yield three clusters of 
the same size with 24 nodes each. Connections within 
the clusters were distributed randomly, and six inter- 
cluster connections were defined to mutually connect all 
clusters. Average density of these networks was again 
similar to the Macaque data. 



2. 2 Methods for detecting important connections 



We tested four candidate measures for predicting vulner- 
able edges in networks. All algorithms were programmed 
in Matlab (Release 12, MathWorks Inc., Natick) as well 
as implemented in C for larger networks. Links were an- 
alyzed as directed connections for all networks. 

First, the product of the degrees (PD) of adjacent 
nodes was calculated for each edge. A high PD indicates 
connections between two hubs which may represent po- 
tentially important network links. 

Second, the absolute difference in the adjacent node 
degrees (DD) of all edges was inspected. A large degree 
difference signifies connections between hubs and more 
sparsely connected network regions which may be im- 
portant for linking central with peripheral regions of a 
network. 



Third, the matching index (MI) ( Hilg etag et al. 2002| ) 
was calculated as the number of matching incoming and 
outgoing connections of the two nodes adjacent to an 
edge, divided by the total number of the nodes' connec- 
tions (excluding direct connections between the nodes 
QSporns 2T )02)). A low MI identifies connections between 
very dissimilar network nodes which might represent im- 
portant 'short cuts' between remote components of the 
network. 

Finally, edge frequency (EF), a measure similar to 
'edge betweenness' IjGirvan and Newman 2002llHolme et al. 20 02). 
indicates how many times a particular edge appears in 
all pairs shortest paths of the network. This measure fo- 
cuses on connections that may have an impact on the 
characteristic path length by their presence in many in- 
dividual shortest paths. We used a modified version of 
Floyd's algorithm IjCormen et al. 2001J) to determine the 
set of all shortest paths and calculate the frequency of 
each edge in it. Multiple shortest paths between nodes i 
and j were present in the analyzed data sets. However, 
the standard algorithm only takes into account the first 
shortest path found. In order to account for edges in al- 
ternative shortest paths, the EF was calculated as the 
average of 50 node permutations in Floyd's algorithm. 
This led to an increased predictive value of this mea- 
sure in all networks; however, the correlations already 
converged after 10 permutations. 

Another possible prediction measure, not used here, 
would be the range of an edge IjWatts 1999l|Sporns~ 2002 ) , 
that is, the length of the shortest path between two adja- 
cent nodes after the edge between them is removed. For 
dense networks, such as cortical connectivity, only range 
values of 2 and 3 occurred. Having only two classes of 
range values was not sufficient to distinguish vulnera- 
ble edges in detail. However, the range may be a useful 
predictor for sparse networks with higher ASP. 
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Table 1 Density, clustering coefficient CC, average shortest path ASP and correlation coefficients r for different vulnerability 
predictors of the analyzed networks (the index refers to the number of nodes) . Tested prediction measures were the product 
of degrees (PD), absolute difference of degrees (DD), matching index (MI), and edge frequency (EF). 





Density 


CC 


ASP 


rpD 


TDD 


TMI 




Macaque73 


0.16 


0.46 


2.2 


0.10** 


0.57** 


-0.40** 


0.84** 


Cat 55 


0.30 


0.55 


1.8 


0.08* 


0.48** 


-0.34** 


0.77** 


AT299 (metabolic) 


0.014 


0.16 


3.5 


0.04 


0.09** 


-0.11** 


0.74** 


EC765 (metabolic) 


0.0067 


0.17 


3.2 


0.31** 


0.38** 


-0.15** 


0.75** 


SC551 (metabolic) 


0.0092 


0.18 


3.3 


0.11** 


0.22** 


-0.04 


0.74** 


SCi846 (protein interactions) 


0.0013 


0.068 


6.8 


0.24** 


0.02 


-0.14** 


0.60** 


German highwayn68 


0.0018 


0.0012 


19.4 


0.19** 


0.06** 


-0.04 


0.63** 


Randorri73 


0.16 


0.16 


1.7 


0.02 


0.06 


0.00 


0.03 


Scale-free73 


0.16 


0.29 


2.0 


0.03 


0.08 


-0.01 


0.03 



* Significant Pearson Correlation, 2-tailed 0.05 level. 
** Significant Pearson Correlation, 2-tailed 0.01 level. 



3 Results 

3. 1 Network patterns underlying edge vulnerability 

The elimination of an edge from a network can have 
two possible effects on the ASP. First, the parts pre- 
viously connected by this edge can still be reached by 
alternative pathways, ff these are longer, the ASP will 
increase. Second, the eliminated edge may be a cut- 
edge, which means that its elimination will fragment the 
network into two disconnected components. The prob- 
ability for fragmentation, naturally, is larger in sparse 
networks. Network separation causes severe damage, as 
interactions between the previously connected parts of 
the network are no longer possible. Therefore, this im- 
pact can be seen as more devastating than the first ef- 
fect, which may only impair the efficiency of network 
interactions. Our ASP calculation disregarded paths be- 
tween disconnected nodes, which would be assigned an 
infinite distance in graph theory. Therefore, the ASP 
in disconnected networks was actually shorter, because 
paths were measured within the smaller separate compo- 
nents. Cut-edges, which lead to network fragmentation, 
frequently occurred in the highway network (30% of all 
edges) and the yeast protein interaction network (23% 
of all edges). However, cut-edges did not occur for cor- 
tical networks of cat and macaque and only to a limited 
extent (<5% of all edges) in the studied metabolic net- 
works. 

In the present calculation both increase and decrease 
of ASP indicate an impairment of the network struc- 
ture, we took the deviation from the ASP of the intact 
network as a measure for damage. We evaluated the cor- 
relation between the size of the prediction measures and 
the damage (Tab^for a ll networks). While most of the 
local measures exhibited good correlation with ASP im- 
pact in real- world networks, the highest correlation was 
consistently reached by the EF measure. For the cortical 
networks, the measures of matching index and difference 
of degrees also show a high correlation. 



Cortical connectivity differed from the other networks 
not only by the performance of different edge vulnerabil- 
ity predictors, but also in the density of connections and 
the amount of clustering. The cortical networks showed 
a higher density than the biochemical metabolic and 
protein-protein interaction networks. Whereas the high- 
way network showed similar density to the biochemical 
networks, its clustering coefficient was much lower be- 
cause of the high proportion of linear paths in the high- 
way network. The random and scale- free benchmark net- 
works — designed to resemble size and edge density of 
the macaque cortical network — are presented at the 
end of the table. 

Fig. ^ shows the ASP after edge elimination plotted 
against edge frequency (EF). For cortical networks (P 
a,b) no network fragmentation occurred, and only an in- 
crease in ASP became apparent. For metabolic networks, 
e.g., the network of 5". cerevisiae jl]c), also a few cut- 
edges, lowering the ASP, were targeted. For the highway 
network containing linear chains of nodes d) , many 
cut-edges were observed. The elimination of these links, 
therefore, resulted in two disconnected compartments, 
each of which had shorter path lengths. 



3.2 Comparison with benchmark networks 

We also calculated the four predictive indices for scale- 
free and random benchmark networks with 73 nodes 
and a similar number of edges as the macaque cortical 
network. Our comparisons focused on this network, be- 
cause it showed the highest correlation between predic- 
tion measures and actual damage for all four measures. 
For the benchmark networks, however, all measures were 
poor predictors of network damage (Fig. [2J . This is sur- 
prising, because scale-free networks generated here by 
growth and preferential attachment appeared to differ in 
their structure from real scale- free networks. Analyzing 
data for one of these real networks, the Internet at the 
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Fig. 2 Evaluation of the performance of four predictors for 
edge vulnerability. Note that the absolute correlation coef- 
ficient was used (MI would have had negative r). For all 
networks, except the random and scale-free benchmark net- 
works, edge frequency had the highest correlation with edge 
vulnerability. In random networks and scale-free networks 
with only one cluster, however, the tested measures were un- 
able to indicate impact of edge elimination. 

autonomous systems level, which was previously shown 
to be scale-free IjBarabasi and Albert 1999|l . we found 
that EF as a prediction measure was also performing 
well in this case (r = 0.62, not shown). The difference 
between the real and simulated scale-free networks may 
result from the fact that scale-free networks generated 
by growth and preferential attachment did not possess 
multiple clusters. We therefore tested whether the lack 
of connections between clusters might be the reason for 
the low performance of EF in the scale-free benchmark 
networks. We generated further 20 test networks; each 
consisting of three randomly wired clusters and six fixed 
inter-cluster connections (Fig. Et) . The inter-cluster con- 
nections (light gray) occurred in many shortest paths 
(Fig. 13b) leading to an assignment of the highest EF 
value, as no alternative paths of the same length were 
available. Furthermore, their elimination resulted in the 
greatest network damage as shown by increased ASP. 



3. 3 Network patterns in biological networks 

After establishing the high impact on ASP of edges with 
large EF, we investigated what made specific edges more 
vulnerable than others. We here discuss two patterns 
that occurred in many of the analyzed networks. First, 
linear chains of nodes that appeared in biological as well 
as the artificial (highway) networks. Second, clusters of 
highly interconnected regions of the network that oc- 
curred for all small- world networks, such as the analyzed 
cortical and biochemical systems. 

Naturally, other and more complex network patterns 
are possible. We already discussed elimination of edges 
between a hub and a node with fewer connections (cf. 
14. ip . Also, the functional role of these patterns (e.g. feed- 
back loops) was not examined here and would merit fur- 
ther study. 
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Fig. 3 Connectivity for multi-clustered benchmark networks 
with comparable density to primate brain connectivity (cf. 
12.1.50 . The gray-level of a connection in the adjacency matrix 
indicates relative frequency of an edge in 20 generated net- 
works. White entries stand for edges absent in all networks, 
a Connectivity of test networks with three clusters and six 
pre-defined inter-cluster connections, b Edge frequencies in 
the all-pairs shortest paths against ASP after elimination of 
edges. Red data points represent the values for the inter- 
cluster connections in all 20 test networks. Inter-cluster con- 
nections not only have the largest edge frequency, but also 
cause most damage after elimination. 



3.3.1 Linear chains Linear chains of nodes with a ter- 
minal end (Fig. 0Ji) became apparent in various biologi- 
cal as well as artificial networks. These patterns were de- 
tected by testing for each node if it was part of a chain, 
that is, if it possessed two edges. In this case the chain 
was followed in both directions, and considered terminal, 
if at least one end of the chain had a terminal node. Using 
this method for identification, each terminal chain con- 
sisted of at least two edges. Nodes in the terminal chain 
were excluded from the further searching process. Ter- 
minal chains occurred frequently in the highway system, 
but also arose in metabolic networks in the form of re- 
dox chain reactions. For the highway system, the average 
terminal chain length was 6.4 edges, with a maximum of 
22 edges. For the yeast protein interaction network, 13% 
of the nodes were part of terminal chains, which were on 
average 2.25 edges long (maximum 6 edges). Eliminat- 
ing edges at the terminal end of a chain would have a 
small impact, as only few nodes become disconnected to 
the rest of the network. On the other hand, severing the 
first edge that connects a chain to the rest of the network 
eliminates all paths leading to the chain nodes from the 
shortest paths matrix. The effect of eliminating edges in 
a chain can be seen clearly for the highway network (Fig. 
nji). Edges that connect chains to the rest of the network 
have a large EF, and their elimination greatly decreases 
ASP, in contrast to edges at the terminal end. Indeed, 
for networks that show many cut-edges, also many ter- 
minal chains occurred. For the highway system, 42% of 
all nodes were part of terminal chains. Similar proper- 
ties of edge vulnerability arise when the terminal of a 
chain end is formed by a small sub-network that is still 
smaller than the main network component at the start 
of the chain. 
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growth that can generate such distributed cluster sys- 
tems ( |Kaiser and Hilgetag 2004| |. 



4 Discussion 




4-1 Measures for identifying vulnerable edges 

We analyzed four measures for identifying vulnerable 
edges and predicting the impact of edge removal on global 
network integrity. Among these measures, the index of 
EF appeared consistently as the best predictor for dam- 
age to edges. The high performance of this measure may 
be linked to characteristic features of biological networks, 
as detailed below Q4.2JI . For the macaque monkey, 7 out 
of the top 10 connections with highest edge frequency 
originated from, or projected to, the amygdala. In addi- 
tion, these edges were most vulnerable as could be ob- 
served by the damage after edge elimination. Therefore, 
the amygdala appears to serve as a central link between 
many clusters of the network. 

Following EF in terms of performance, the index for 
the difference of degrees also showed a high correlation in 
both the cortical and metabolic networks (Fig. [2J . This 
means that connections between highly and sparsely con- 
nected nodes are vulnerable, especially in cortical net- 
works. When a node with few connections is connected 
3.3.2 Clustered architecture A clustered or modular ar- to an already well-connected node (hub), it can access 
chitecture is a characteristic feature of many naturally a large part of the network, by using routes involving 
occurring networks, such as cortical connectivity net- the hub. After eliminating the connecting edge, the node 
works in the pri mate fVoung 1992||Young 1993||Hilgetag et aiyc3IDI@0) iave to use longer alternative pathways to reach 
or the cat brain ljScannell et al. 1995IIScannell et al. 19991 the same parts of the network. This effect is particularly 
Hilg etag et al. 2000| ) as well as metabolic networks QRavasz etete ioi2fl,02).the node was the only one in its local neigh- 
These systems are known to consist of several distinct, borhood that was connected to the hub. 



Fig. 4 Network pattern and corresponding edge vulnerabil- 
ity, a Elimination of edges of linear chains of nodes results 
in two disconnected components and a lower ASP. Edges 
eliminated at the proximal end of the chain (red) cause a 
larger change in ASP than at the terminal end (blue), b 
For clustered networks, edges within the clusters (blue) can 
be replaced by several alternative pathways. Therefore, their 
elimination causes a smaller increase of ASP than that of 
edges between clusters (red). 



linked clusters with a higher frequency of connection 
within than between the clusters. Inter-cluster connec- 
tions have also been considered important in the context 
of social contact networks, as 'weak ties' between indi- 
viduals ( Granovet ter 1973|l and separators of communi- 
ties l|Girvan and Newman 2 002). We, therefore, specu- 
lated that connections between clusters might be gen- 
erally important for predicting vulnerability (Fig. 0Jd). 
Whereas many alternative pathways exist for edges within 
clusters, the alternative pathways for edges between clus- 
ters may be considerably longer. Interestingly, previously 
suggested growth mechanisms for scale-free networks, 



The matching index showed a large (negative) corre- 
lation with edge vulnerability in cortical networks and — 
to a lower degree — the protein-protein interaction net- 
work. Therefore, edges between dissimilar connections, 
that is, with low MI, were more vulnerable. This was due 
to the cluster structure of these networks (cf . 14.2(1 . Nodes 
with similar connectivity belong to the same cluster, 
and therefore multiple alternative pathways are avail- 
able. Dissimilar nodes are more likely to be part of dif- 
ferent clusters with few alternative pathways. 

From theoretical studies it has been proposed for 
scale-free networks that edges between hubs are most 



such as preferential attachment (Barabas i and Albert 1999jl . vulnerable ((Holme et al. 200*2")) . However, in the networks 
or strategies for generating hierarchical networks (Barabasi ctaai]al20Qd ) here, both in the scale- free yeast protein in- 



did not produce distributed, interlinked clusters. Con- 
sequently, the low predictive value of EF in the scale- 
free benchmark networks was attributable to the fact 
that scale-free networks grown by preferential attach- 
ment consisted of one central cluster, but did not possess 
a multi-cluster organization. This suggests that alterna- 
tive developmental models may be required to reproduce 
the specific organization of biological networks, and we 
have recently presented an algorithm based on spatial 



teraction network as well as for cortical networks, edges 
which connected nodes possessing many connections (large 
product of degrees) were not particularly vulnerable. Al- 
though a low increase can be seen for the correlation co- 
efficient, edges with maximum vulnerability occurred for 
a small product of degrees. 

One of the general advantages of using prediction 
measures, instead of testing the damage for all edges of a 
network, is computational efficiency. For a network with 
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e edges and n nodes, the order of time for the calcula- 
tion of EF that can predict the effect of edge elimination 
for all edges is 0(n 3 ). This is lower than for testing the 
damage after edge elimination for all edges individually, 
calculating the ASP e times resulting in a time complex- 
ity of 0(e • n 3 ). The prediction measures presented here 
might therefore be particularly useful for large networks 
in which global testing is computationally impractical, 
or for networks with frequently changing connections de- 
manding a regular update. Examples of these include ac- 
quaintance networks, Internet router tables, and traffic 
networks. Identifying vulnerable edges by EF generally 
appears to work well for various biological networks. 

4-2 Vulnerable edges in biological networks 

In the analyzed biological networks, inter-cluster projec- 
tions may play an important role in linking functional 
units. For cortico-cortical networks, they connect and 
integrate different sensory modalities (e.g., visual, audi- 
tory) or functional sub-components ( jHilgetag et al. 2000| ) 
A lesion affecting these connections may result in disso- 
ciation disorders IjGeschwind 1965|l . 

For metabolic networks, reactions proceed more fre- 
quently within a reaction compartment (e.g., mitochon- 
dria and endoplasmatic reticulum) than between com- 
partments. Therefore, in these systems as well, local- 
ized clusters arise, with many reactions within a com- 
partment and few connections between compartments. 
Such an organization is also found in the investigation of 
protein-protein interactions and their spatial and func- 
tional clustering, in which fewer proteins from different 
groups interact IjSchwikowski et al. 2000(1 . Once again, 
interactions between proteins from different compart- 
ments correspond to inter-cluster connections, and might 
thus be among the most vulnerable edges of the network. 

To a lower extent, 'inter-cluster' connections also ap- 
peared in the highway network, as the highway sub- 
systems of Western and Eastern Germany formed (spa- 
tially separate) dense regions connected by only four 
highways. That is, the elimination of four edges would 
once again split the German highway system into East- 
ern and Western components. Our analysis of artificial 
networks was restricted to networks without functional 
differentiation of edges or nodes. It remains to be seen if 
functional or social networks, for instance, interactions 
of people with different functions within a company, may 
show a higher similarity in cluster-architecture with bi- 
ological networks. 

Clustered network architecture appears to result in 
edge-robustness in a similar way as scale-free architec- 
ture results in node-robustness. Random elimination of 
edges will most frequently select edges within a cluster. 
For paths routed through these edges, various alternative 
pathways exist, and the damage after edge elimination 
is small. Targeted attacks on inter-cluster connections, 



on the other hand, result in large network damage. We 
note that the cortical networks investigated here exhibit 
properties of both small- world ( Hilgc tag et al. 2000| ) and 
scale- free (IKaiser et al. 2004|) networks and are therefore 
particularly robust to random failure of edges and nodes. 

4-3 Conclusions 

We examined how local network features and patterns 
relate to global network properties, such as robustness 
towards edge elimination. The correlation between dif- 
ferent predictors and the actual damage after the elimi- 
nation of a connection signified differences in the global 
network architecture. Specifically, various biological net- 
works appear to be organized as distributed, linked net- 
work clusters. We have shown that inter-cluster connec- 
tions represent the most vulnerable edges in these net- 
works, and that their position can be predicted using the 
edge frequency measure. 
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